Nonlocal Andreev reflection at high transmissions 
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We analyze non-local effects in electron transport across three-terminal normal-superconducting- 
normal (NSN) structures. Subgap electrons entering S-electrode from one N-metal may form Cooper 
pairs with their counterparts penetrating from another N-metal. This phenomenon of crossed An- 
dreev reflection - combined with normal scattering at SN interfaces - yields two different con- 
tributions to non-local conductance which we evaluate non-perturbatively at arbitrary interface 
transmissions. Both these contributions reach their maximum values at fully transmitting interfaces 
and demonstrate interesting features which can be tested in future experiments. 



At sufficiently low temperatures Andreev reflection 
(AR) [1] dominates charge transfer through an inter- 
face between a normal metal and a superconductor (NS): 
An electron propagating from the normal metal with en- 
ergy below the superconducting gap A enters the super- 
conductor at a length of order of the superconducting 
coherence length ^, forms a Cooper pair together with 
another electron, while a hole goes back into the nor- 
mal metal. As a result, the net charge 2e is transferred 
through the NS interface which acquires non-zero subgap 
conductance 
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In hybrid NSN structures with two N-terminals, elec- 
trons may penetrate into a superconductor through both 
NS interfaces. Provided the superconductor size (dis- 
tance between two NS interfaces) L strongly exceeds ^, 
AR processes at these interfaces are independent. If, 
however, the distance L is smaller than or comparable 
with ^, two additional non-local processes come into play 
(see Fig. 1). Firstly, an electron with subgap energy 
propagating from one N-metal can penetrate through the 
superconductor into another N-electrode with the prob- 
ability ^ exp(— L/^). Secondly, an electron penetrating 
into the superconductor from the first N-terminal may 
form a Cooper pair by "picking up" another electron from 
the second N-terminal. In this case a hole will go into 
the second (not the first!) N-metal and, hence, AR turns 
into a non-local effect. The probability of this process 

- usually called crossed Andreev reflection (CAR) [J 01 

- also decays as ~ exp(— L/^) and, in combination with 
direct electron transfer between normal electrodes, de- 
termines non-local conductance in hybrid multi-terminal 
structures which can be directly measured in experiment. 

CAR has recently become a subject of intensive in- 
vestigations both in experiment [^, I^, 0] and in theory 
d, [9I lid. 11 , III (see also further references therein) . Al- 



FIG. 1: Two elementary processes contributing to non-local 
conductance of an NSN device: (1) direct electron transfer 
and (2) crossed Andreev reflection. 



perconducting critical temperature [1,|3|) as well as zero- 
bias anomalies in the Andreev conductance due to both 
disorder-enhanced interference of electrons 
and Coulomb effects 
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though a non-local conductance was observed in all these 
experiments, an unambiguous and detailed interpretation 
of the existing experimental data still remains a chal- 
lenge, to a certain extent because in addition to the above 
processes a number of other physical effects may consid- 
erably influence the observations. Among such effects we 
mention, e.g, charge imbalance (relevant close to the su- 



CAR is also sensitive 
to magnetic properties of normal electrodes. Although 
theoretical investigation of the above physical effects is 
certainly of interest and may help to account for some ex- 
perimental observations, we believe that, beforehand, it 
is important to reach quantitative understanding of CAR 
in simpler situations when (at least some of) the above 
effects can be disregarded. 

As in most cases metallic interfaces are not fully trans- 
parent, AR is usually combined with normal electron 
scattering at such interfaces. The relative "weights" of 
these two processes are determined by interface trans- 
mission. In the case of multi-terminal hybrid structures 
normal reflection, tunneling, local AR and CAR combine 
in a complicated and non-trivial manner. For instance, 
it was demonstrated [1, Q that in the lowest order in the 
interface barrier transmission and at T = CAR con- 
tribution to cross-terminal conductance is exactly can- 
celled by that from elastic electron cotunneling [3] , while 
no such cancellation is expected in higher orders in the 
transmission [10]. However, complete theory of non-local 
phenomena in question which would fully describe an in- 
terplay between all scattering processes to all orders in 
the interface transmissions and set the maximum scale 
of the effect remains unavailable. Such a theory requires 
non-perturbative methods and is the main subject of the 
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FIG. 2: Schematics of our NSN device. 
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present work. 

The model and formalism. Consider three-terminal 
NSN structure depicted in Fig. 2. We will assume that all 
three metallic electrodes are non-magnetic and ballistic, 
i.e. the electron elastic mean free path is large. Transmis- 
sions Di and D2 of two SN interfaces (with cross-sections 
Ai and A2) may take any value from zero to one. The 
distance between the two interfaces L as well as other ge- 
ometric parameters are assumed to be much larger than 
y^Ai,2, i.e. effectively both contacts are metallic con- 
strictions. In this case the voltage drops only across SN 
interfaces and not inside large metallic electrodes. Hence, 
nonequilibrium (e.g. charge imbalance) effects related to 
the electric field penetration into the S-electrode can be 
neglected. In what follows we will also ignore Coulomb 
effects [iiHIiTl- 

For convenience, we will set the electric potential of 
the S-electrode equal to zero, V = 0. In the presence of 
bias voltages Vi and V2 applied to two normal electrodes 
(see Fig. 2) the currents Ii and I2 will flow through SNi 
and SN2 interfaces. These currents can be evaluated with 
the aid of the quasiclassical formalism of noneq uilibrium 
Green-Eilenberger-Keldysh functions g^'^'^ [ig|. In the 
ballistic limit the corresponding equations take the form 



£f3 + eV{r, t) - A(r, t), g"^'^'"" {pp, e, r, t) 



^^;FV.9^■'^■^(pF,e,r,^) = 0, 



(1) 



where [a, b] = ah — ba, e is the quasiparticle energy, Pf = 
mvp is the electron Fermi momentum vector and fa is 
the Pauli matrix. The functions g^'^'^ also obey the 
normalization conditions (5^)^ = (ff"^)^ = 1 and g^g^ + 
gK^A _ g jjgj-g below the product of matrices is 
defined as time convolution. 

The matrices q and A have the standard form 
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R,A,K 



gR,A,K jB,A,K 
jR,A,K gB.,A.K ] ' 



A = 



A 

-A* 



(2) 



where A is the BCS order parameter. The current den- 
sity is related to the Keldysh function g^ as 

jir.t) ^ -'^ j de{vFS^[h9^{PF,e,r,t)]) , (3) 



FIG. 3: Quasiclassical trajectories contributing to local (a) 
and non-local (b and c) currents. 



where Nq = mpF/^Tr'^ is the density of state at the Fermi 
level and angular brackets (...) denote averaging over the 
Fermi momentum directions. 

The above equations should be supplemented by ap- 
propriate boundary conditions. In order to match quasi- 
classical Green functions at the N- and S-sides of SNi in- 
terface (respectively gNi and gs) we will make use of Za- 
itsev boundary conditions [20| for matrices 9 — (^^q ^-a 



9" 



■9n, 



9ni = .9s 



3s ' 



r[Ri{tr + i9-r] = D^g-g^ 



(4) 
(5) 



where g^ = + ±g^ ±gs, g^ = gs{±Px) (see 
Fig. 3a), Ri{pxi) = I- Di{px^),Pxi is the component of 
Pf normal to the SNi interface. Green functions at SN2 
interface are matched analogously. Deep inside metallic 
electrodes S, Ni and N2 the Green functions should ap- 
proach their equilibrium values g^'^ = ±(£73 — 
in a superconductor and g^^^ = ifs in normal metals, 
n^'^ ~ ^{e ± iSy — A^. For the Keldysh functions far 

from interfaces we have g^ — 5^ ( V ft"- ) ^ ( V h- ) 9^^ 
where h± = tanh[(e±ey)/2T]. Voltage in above expres- 
sion equals to F = 0, Vi and V2 respectively in S, Ni and 
N2 electrodes. The parameter A is chosen to be real. 

Relevant trajectories. Electron trajectories which con- 
tribute to the current Ii through SNi interface are shown 
in Fig. 3. Trajectories presented in Fig. 3a do not enter 
the terminal N2 and yield the standard BTK contribu- 
tion to /i. In addition there exist trajectories (Fig. 
3b, c) involving all three electrodes. They fully account 
for all scattering processes - both normal and AR - to all 
orders in the interface transmissions and determine non- 
local conductance of our NSN device. As follows from 
Fig. 3b, c for each direction of px one can distinguish 
four different contributions to non-local conductance cor- 
responding to different trajectory combinations. 

Note that applicability of the above quasiclassical for- 
malism with boundary conditions ([5]) to hybrid structures 
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with two (or more) barriers is, in general, a non-trivial 
issue 21[ which requires a comment. Electrons scattered 



at different barriers may interfere and form bound states 
(resonances) which cannot be correctly described within 
a formalism employing Zaitsev boundary conditions [20| . 
In our geometry, however, any relevant trajectory reaches 
each interface only once whereas the probability of multi- 
ple reflections at both interfaces is small in the parameter 
AiA2/L'^ <C 1. Hence, resonances formed by multiply 
reflected electron waves can be neglected, and our for- 
malism remains adequate for the problem in question. 

Quasiclassical Green functions. The above equations 
can be conveniently solved introducing parameterization 
of the matrix Green functions g^-^-^ by four Riccati am- 
plitudes and two "distribution functions" [121 ■ This pa- 
rameterization allows to transform Eq. ([1]) to a set of 
decoupled equations. It is also important that non-linear 
Zaitsev boundary conditions ([?]), ([5]) can be rewritten 
in terms of Riccati amplitudes and "distribution func- 
tions" in a rather simple form [22]. Integration of the 
resulting equations along the trajectories shown in Fig. 
3 is straightforward. Finally we arrive at the following 
expression for the Keldysh Green function g^_^ at SNi 
interface (on the N-metal side) 



(6) 



Here g^ai^i) comes from the trajectories of Fig. 3a re- 
sponsible for the BTK current at SNi interface, while two 
other terms come from the trajectories of Fig. 3b, c which 
also involve N2-electrode. The term 5ff,_|_j,(V'i) yields a 
correction to the BTK term which will be discussed later. 
The last contribution h^^i^^) accounts for non-local 
conductance of our device. For positive p^^ > we have 



9l2,b+ciV2) = 2DiD2 



1 — tanh^ iLQ/vp 
P{Ri,R2) 



2T 



6'fci?2|aptanh 



2T 



+ 9cRi\a\^ tanh ^ + 9b tanh ^ ] , (7) 



2T 



2T 



where we defined fl = Q^, P(i?i,i?2) = |1 - i?ii?2a^ - 
Q[e{l+RiR2a^)+Aa{Ri+R2)]\^, Q = tanhim/vF, 
a — (fl — e)/A, 9b and 9c equal to unity for trajectories 
of respectively Fig. 3b and 3c and to zero otherwise. As 
expected, Eq. ([7]) identifies four different contributions 
entering with the corresponding amplitudes and reflec- 
tion coefficients. Note that only one out of these con- 
tributions survives in the case of rcflectionless interfaces. 
In contrast, for weakly transmitting barriers (i?i,2 ~^ 1) 
and e < A all four terms enter with equal prefactors. 

As for the function , at SNi interface it does not 
depend on V2 for positive p^^ > 0. The values of g^ and 
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FIG. 4: (Color online) Differential non-local conductance at 
T = as a function of voltage for Di — 0.5, D2 = 0.8 and 
different L. Inset: the same for eV < A. 
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g" for negative Pxi < are easily recovered by means of 
the relation g^{—pF, — £, r, t) = g^{pF, £, r, t). 

Non-local conductance. Substituting the results ©, 
into Eq. (O we obtain 



/i -/ll(^^l) + /l2(T4), 
/2 -/2l(^^l)+/22(V^2). 



(8) 
(9) 



Here In and /22 consist of the standard BTK currents 
0, [lOl and CAR terms to be specified later and 



h2{V)=l2l{V) = - 



G 



N12 



2e 



de 



, e + eV , e ' 

tanh tanh — 

2T 2T 



^ , ,9^,. ^ , ,,,1 — tanh iLfl/vp , „x 
X (l-7^l|a^(l-7^2|a^ ^^^^ ^ , (10) 

where Vy,2 = 1 - '^1,2 = R>\,2{VFlia) and pf1i(2) is 
normal to the first (second) interface component of the 
Fermi momentum for electrons propagating straight be- 
tween the interfaces, 



G 



N12 



87i72A/"iA/"2PiP2 

RqPlL^ 



(11) 



is the non-local conductance in the normal state, J\fi.2 = 
PF-^i,2/47r define the number of conducting channels of 
the corresponding interface, Rq — 27r/e^ is the quantum 
resistance unit. Eq. I|10p represents the central result 
of our paper. This expression fully determines non-local 
conductances of our NSN device at arbitrary transmis- 
sions of SN interfaces. 

The differential non-local conductance evaluated with 
the aid of Eq. (fTO|) at T = is presented in Fig. 4 at 
sufficiently high interface transmissions. We observe that 
this quantity increases sharply around eV A and ap- 
proaches the L-independent (normal) limit at eV S> A. 
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In the limit T, Vi^2 ^ A only subgap quasiparticles con- 
tribute and the differential conductance becomes voltage- 
independent. We have /12 = — G12V2, where 



G12 _ ViV2il-tanh^ LA/vf) 

Gni2 ~ [1 + 7^l7^2 + {Til + 7^2) tanh LA/vp] 



(12) 



The value G12 (HH) gets strongly suppressed with de- 
creasing X'i_2 and increasing i, as also seen in Fig. 4. 
Note, that the dependence of G12 on L reduces to purely 
exponential at all L only in the lowest nonvanishing or- 
der in the transmission of at least one of the barriers, 
e.g., G12 oc D^P^ exp(-2iA/wF) for Di,2 < 1, whereas 
in general this dependence is slower than exponential 
at smaller L and approaches the latter only at large 
L > wf/A. 

For a given L the non-local conductance reaches its 
maximum in the case of reflectionless interfaces Di 2 = 1- 
Interestingly, in this case for small L <C vp/A the con- 
ductance G12 identically coincides with its normal state 
value Gni2 a-t any temperature and voltage. This re- 
sult can easily be understood bearing in mind that for 
£'1,2 = 1 only trajectories indicated by horizontal lines 
in Fig. 3b, c contribute to G12. For i — » there is 
"no space" for CAR to develop on these trajectories and, 
hence, CAR contribution to G12 vanishes, whereas di- 
rect transfer of electrons between Ni and N2 remains 
unaffected by superconductivity in this limit. 

The situation changes provided at least one of the 
transmissions is smaller than one. In this case scattering 
at SN interfaces mixes up trajectories connecting Ni and 
N2 terminals with ones going deep into and coming from 
the superconductor. As a result, CAR contribution to 
G12 does not vanish even in the limit L and G12 
turns out to be smaller than Gni2- 

Finally, we would like to briefly address the non-local 
correction to Gn which arises from the CAR process de- 
scribed by the term ^^^^^(Fi) in Eq. dH). At T, Vx^i < A 
we have /n — GuVi, where Gn = Gf + ^G\\. Here 
Gf^^ is the standard BTK term 



qBTK ^ 

Rn 



(13) 



and for the non-local term we obtain 



2(1 + 7^2)(l - tanh^ LA/vp) 



SGn ^ 

Gn^2 [1 + 7^l7^2 + (Tei +7^2)tanhLA/l;F]2 
Vi [(1 + 7^2 tanhLA/uF)^ + 3(7^2 + tanhLA/vpf] 



V2[l + 7^l7^2 + (7^l + 112) tanhLA/uF]^ 



(14) 



As compared to the BTK conductance (|13p the CAR 
correction ([T^ contains an extra small factor A2/L'^ 
and, hence, in many cases can be neglected. On the 
other hand, since CAR involves tunneling of one electron 



through each interface, for small /?i <C 1 and D2 « 1 we 
have 5Gii cx i.e. for Di < {A2/ L'^) exp{-2LA/vp) 
the CAR contribution (fT4|l may well exceed the BTK 
term Gf cx 

In summary, we have developed a theory of non-local 
electron transport in ballistic NSN structures with arbi- 
trary interface transmissions. Non-trivial interplay be- 
tween normal scattering, local and non-local Andreev re- 
flection at SN interfaces yields a number of interesting 
properties of non-local conductance which can be tested 
in future experiments. 

We would like to thank V. Chandrasekhar for commu- 
nicating the results to us prior to publication and to 
A. A. Golubov for useful discussions at an early stage of 
this work. 
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